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The Kane-Mele (KM) model is proposed to describe 
the quantum spin Hall effect (QSHE) of electrons on 
the two-dimensional honeycomb lattice.^' ^ Here, we will 
show that, in a certain parameter region, the London 
equation is obtained from the effective field theory of 
the KM model with an electronic correlation.'^"^ Wc use 
h = c = 1 unit and the Minkovskian metric g^'' = 
diag{l, —1, —1), where /i, z^ = 0,x, y. 

We consider a layered honeycomb lattice system in 
which each layer is described by the KM model. We as- 
sume that interlayer coupling is negligibly small. One of 
the essential ingredients of the KM model^'^ is intrinsic 
spin-orbit (SO) coupling Aso, which generates an effec- 
tive magnetic field depending on spin as well as an ex- 
citation gap A = 3-\/3Aso to the electron.^ Thus, the 
term enables quantization of the spin Hall conductivity 
(SHC),i'2 
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where d is the interlayer distance. The model can also 
have the Rashba extrinsic SO coupling A^, which breaks 
the inversion symmetry and is induced by an electric 
field perpendicular to the honeycomb lattice plane. The 
term also breaks the conservation of electron spin Sz/2. 
Hereafter, we consider the case of Xr = 0. 

We add the on-site Coulomb repulsion U > Q. The 
Hamiltonian per a layer is 
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where c; {c]) is the annihilation (creation) operator of 
an electron with spin at the i-th site and t is nearest 
neighbor hopping. The second term is the intrinsic SO 
term consisting of next nearest neighbor hopping, and 
Vij = -j={di X d2)z = il, where di and ^2 are unit 
vectors along the two bonds where the electron moving 
from site j to i passes. 

Let us discuss how to deal with the electron correlation 
U. On-site coulomb repulsion can be written by the on- 
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site spin-spin interaction: 

Uni^n,i = ^('M + "-4) ~ T-(4*Cj)^- (3) 

The first term merely gives the renormalization for the 
chemical potential and can be neglected. We introduce 
the auxiliary field (pi, which is a three-component vec- 
tor in the spin space, and use the Stratonovich-Hubbard 
transformation,^ H — > Hsh = H + AH, where 
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The spin-spin interaction is eliminated in appearance, 
but we have coupling between Lpi and the electron spin, 
and a quadratic term of Lpi instead. 

We consider the continuum limit and take into account 
the low-energy electronic excitations around K and K' 
points in the Brillouin Zone.^'^ We introduce the elec- 
tromagnetic C/(l) gauge field A^ and SU{2) spin gauge 
field d^ via the covariant derivative. 
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where do = ip (the auxiliary field in the continuum limit) 
and a is a constant external field introduced artificially 
to estimate the spin current. We define a parameter 

AUa'^d 
5-^^, (7) 



where a is the lattice constant, and the microscopic La- 
grangian density is 

C ^ ¥ {iDo - iv{D^ 

2 2^lo 

"^Tcrs is the eight-component Fcrmion field 
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where ^ 

labeled by the eigenvalues of the diagonal components of 
valley spin t. sublattice spin a and real spin s/2. The 
parameter v is the Fermi velocity when the system is in 
the metallic state, and eo and ^0 denote the dielectric 
constant and magnetic permeability, respectively. Note 
that, except for the last term, the Lagrangian (8) pos- 
sesses the C/(l)cm X U{l)z local gauge symmetry. The 
SU{2) gauge symmetry is broken down to U{l)z, since 
the SO term contains s^. 

The calculation shown below is similar to that pre- 
sented in Ref. 7, although the physical meaning of the 
spin gauge field is different. We integrate out Vf and 
obtain the one-loop effective Lagrangian for the gauge 
fields in the low-energy (long-wavelength) region com- 
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The result is^ 
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+ (tcrnis independent of a^^ and A^J, (10) 

where £ind stands for the induced part of the effective 
Lagrangian. The first term in Eq. (10) is the BF term.^ 
The coefficient is the quantized SHC given in eq. (1). 
Note that only a^ couples to the electromagnetic gauge 
fields, because the SU{2) symmetry is broken down to 
the U{l)z symmetry by SO coupling. The Maxwell term 
is renormalized as'''^ 
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By using the parameters in Table I and the relations 
■^ ~ 1/137 and eoA^o = 1, we obtain eo/ g^Aid "= ^-^ 

2 2 o 

and ^iq ■ g^i^i^ = 2 X 10" , i.e., /i ~ /io- The elastic term 
for «§(= <p^) is also induced. We can recognize that any 
potential terms (i.e., zeroth-order terms with respect to 
the derivative 9^) of A^ and also a^ in £ind are absent 
because of the presence of C/(l)om x U{l)z gauge symme- 
try in the Fermionic part of the microscopic Lagrangian 
(8). Thus, the low-energy and long- wavelength physics of 
Afj^ and a^ is described definitely by eq. (9). 

We consider the static magnetic response. The equa- 
tions of motion obtained from eq. (9) are 
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If we obtain B from eqs. (13) and (14), we can find that 
the term ag/(? in eq. (13) can be neglected when 
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Going back to the microscopic Lagrangian (8), we may 
disregard the quadratic term in this case, and thus, a§ 
can be regarded as "the spin chemical potential" . 

Here, we make two crucial assumptions: (i) eq. (15) is 
satisfied and (ii) the system is not in the topological Mott 
insulating phase^ but in the topological band insulating 
phase, which exhibits the quantum spin Hall effect and 
its low-energy and long-wavelength physics is described 
by the effective Lagrangian eq. (9) when the spin is con- 
served (e.g., Afl = 0). To realize (ii), U should not be too 
large. ^'^ Thus, a large a is favorable to realize (i) and 
(ii) simultaneously [see eq. (7)]. Note that a large d is 
unfavorable since a^ is suppressed [see eq. (1)]. 

As mentioned, the 1/g term in eq. (13) can be ne- 
glected under condition eq. (15). The r.h.s. of eq. (13) is 
a result of the fact that the spin density is proportional 
to the magnetic field B owing to the BF term. It resem- 
bles the Zeeman effect, but the essential difference is that 
the coefficient is not the Bohr magneton but the SHC. 
The r.h.s of eq. (14) is a result of the fact that, owing 



to the BF term, the electric current fiows perpendicular 
to the gradient of the spin chemical potential Og. This 
may be called the dual quantized spin Hall effect (dual 
QSHE). Taking the rotation of both sides of eq. (14) and 
using eq. (13) (note, again, that we are neglecting the 
1/g term), we obtain the London equation,^' ^*^ 
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Let us check condition (15). By using the parameters for 
the topological band insulator Na2lr03'^ in Table I, we 
obtain crj^/ig — 7.0 x 10"^, i.e., eq. (15) is not satisfied. 
We note that, if we have a superlattice structure with 
a = lO^'A, eq. (15) is satisfied since cr^^^g — 7. 

The physical implication of eq. (16) is the Meissner 
effect. The penetration depth of the magnetic field is 
estimated using the parameters in Table I as Apen. = 
27r(2ecrJj^)"Mei/2/^-i/2 ^ 3990 1, which is as short as 

the typical value of the order of 1000 A for superconduc- 
tors. It should be pointed out that the large A ~ 0.5 eV 
enhances the Meissner effect, since ^e ex 1/A. 

To discuss the Meissner effect more precisely, the sam- 
ple boundary should be taken into account. A further 
discussion of the role of the helical edge states^' ^ in elec- 
tron correlation is needed. 



Table I. Parameters used for estimations. These are typical 
values for Na2lr03,^ which is a honeycomb-layered topological 
band insulator with electron correlation. 
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U 
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V 


0.5cV 


0.5eV 


lOA 


lOA 


3 X 10*m/s 
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